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ABSTRACT 

Aims. The luminous blue variable (LBV) stars are peculiar very massive stars. The study of these 
stellar objects and their surroundings is important for understanding the evolution of massive 
stars and its effects on the interstellar medium. We study the LBV star candidate G26.47+0.02. 
Methods. Using several large-scale surveys in different frequencies we performed a multiwave- 
length study of G26.47+0.02 and its surroundings. 

Results. We found a molecular shell (seen in the 13 CO J= 1—0 line) that partially surrounds the 
mid-infrared nebula of G26.47+0.02, which suggests an interaction between the strong stellar 
winds and the molecular gas. From the HI absorption and the molecular gas study we conclude 
that G26.47+0.02 is located at a distance of ~ 4.8 kpc. The radio continuum analysis shows a both 
thermal and non-thermal emission toward this LBV candidate, pointing to wind-wind collision 
shocks from a binary system. This hypothesis is supported by a search of near-IR sources and 
the Chandra X-ray analysis. Additional multiwavelength and long-term observations are needed 
to detect some possible variable behavior, and if that is found, to confirm the binary nature of the 
system. 

Key words. Circumstellar matter - stars: massive - stars: mass-loss - stars: individual: 
G26.47+0.02 - ISM: clouds 

1. Introduction 

Luminous blue variable (LBV) stars are peculiar very massive stars that evolve from the O-type 
main sequence burning hydrogen in their core to become a Wolf-Rayet (WR) helium core burning 
star. Their main characteristics are a high mass-loss rate (up to 10 -4 M yr _1 ), sometimes accompa- 
nied by so-called giant eruptions, a high luminosity (~ 10 6 L Q ), and significant photometric as well 
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as spectroscopic variability (Humphreys & Davidson 1994). The high mass-loss rate associated 
with the LBV phase typically results in the formation of an ejecta nebula around the star (e.g. Nota 
et al. 1995; Clark et al. 2003). The nebulae around LBV stars generally are strong emitters in the 
mid-infrared, showing that they are composed of both gaseous and dusty components. 

At present, there are very few cases in which molecular material related to LBV or LBV- 
candidate nebulae have been determined by molecular studies. The most representative sources are 
AG Car (Nota et al. 2002) and G79.29+0.46 (Rizzo et al. 2008; Jimenez-Esteban et al. 2010). Very 
recently, Petriella et al. (2012) discovered a fragmented molecular shell delineating the infrared 
bipolar outer shell of the LBV star G24.73+0.69. The authors argued that the molecular shell 
formed from the interstellar material that was swept-up by the stellar wind of the central star. 
On the other hand, they detected some molecular emission probably associated with the inner 
G24.73+0.69 infrared nebula, which may have originated in a mass-ejection event from the central 
star. 




Galactic Longitude 



Fig. 1. Color-composite image of the G26 surroundings, where the Spitzer-IRAC 8 pm emission is displayed 
in green, the Spitzer-MIPSGAL emission at 24 yum in red, and the radio continuum emission at 20 cm is 
presented in blue with white contours with levels of 2.1, 4.1, and 6.5 mjy beam -1 . 

G26.47+0.02 (hereafter G26) is an LBV star candidate. This source presents a compact 
mid-infrared nebula, which closely resembles the ring nebulae around the LBV star candidates 
G79.49+0.26 and Wra 17-96 (Wachter et al. 2010; Clark et al. 2003; Egan et al. 2002). According 
to Clark et al. (2003), adopting a distance of 6.5 kpc, the G26 central star has a luminosity of 10 6 
L and a mass-loss rate of 9 x 10~ 5 M yr _1 , which indicated that it could be one of the most ex- 
treme stars in the Galaxy, similar to the known LBVs AG Car and AFGL 2298. Additionally, the 
authors suggested that G26 is photometrically variable. Very recently Naze et al. (2012) performed 
an X-ray survey of Galactic LBV stars. The authors reported the detection of G26 by Chandra and 
concluded that the X-ray emission is well explained by wind- wind collisions in a possible binary 
system. 

In Figure 1 we present the LBV star candidate G26 in a color-composite image with the Spitzer- 
IRAC 8 pm emission (green), Spitze r-MIPSGAL emission at 24 pm (red), and the radio continuum 
emission at 20 cm (blue with white contours). We used the mosaicked image from GLIMPSE in 
the Spitzer-IRAC band at 8 pm, which has an angular resolution of ~179 (see Fazio et al. 2004 
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and Werner et al. 2004). MIPSGAL is a survey of the same region as GLIMPSE, using the MIPS 
instrument (24 and 70 /urn) on Spitzer. The MIPSGAL resolution at 24 jum is 6". On the other hand, 
the radio continuum data at 20 cm with a FWHM synthesized beam of about 5" was extracted from 
the New GPS of the Multi-Array Galactic Plane Imaging Survey (Helfand et al. 2006). It can be 
appreciated from Fig. 1 that the nebula surrounding the central star brightens at 24 fim and presents 
weak emission at 8 /mi, a different case compared with the LBV star G24.73+0.69, which is very 
bright in both bands. In addition, G26 has associated radio continuum emission, as was previously 
reported by Clark et al. (2003). 

Bellow we study the interstellar medium around G26 and perform a multiwavelength analysis 
toward the central stellar object to unveil its nature. 

2. The environment of the LBV star candidate G26 

2.1. HI absorption: distance estimate 

In this section we use HI data extracted from the VLA Galactic Plane Survey (VGPS; Stil et al. 
2006), which have an angular and spectral resolution of ~V and 1.56 km s _1 , respectively. 




Fig. 2. HI spectra obtained toward G26.47+0.02 and its surroundings. The spectrum obtained toward the 
source (the On position) is presented in blue, in red is presented the averaged HI emission taken from four 
positions separated by approximately two beams from the source in the direction of the four galactic cardinal 
points (the Off position), and the subtractions between them is presented in black. The 3cr uncertainty of the 
subtraction is ~8 K, which is shown with the dashed line. 

Taking into account that G26 presents radio continuum emission at 20 cm, it is possible to 
estimate its distance through a study of the HI absorption and the related molecular gas. Figure 
2 shows the HI spectra toward the source and surroundings. The HI emission obtained over the 
source (the On position: a beam over the radio maximum of the source) is presented in blue, in 
red is presented the average HI emission taken from four positions separated by approximately 
two beams from the source in direction of the four galactic cardinal points ( the Off position), 
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and the subtraction between them is presented in black, which has a 3<x uncertainty of ~8 K. 
The figure shows that the last absorption feature appears at v L sr ~ 66 km s _1 . From the Galactic 
rotational model of Fich et al. (1989) this velocity implies the kinematic distances of either 4.3 kpc 
or 11.0 kpc. Taking into account that the tangent point (at Vlsr ~ 119 km s _1 ) does not present 
any absorption, following Kolpak et al. (2003), we favor the near kinematic distance. Considering 
that Clark et al. (2003) used 6.5 kpc as a maximum distance to the source, we conclude that G26 
should be located between 4.3 and 6.5 kpc. In the following section, based on this result and by 
considering the presence of related molecular gas, the distance is constrained with better precision. 

2.2. Molecular gas 

We use molecular data extracted from the Galactic Ring Survey (GRS). The GRS was performed 
by the Boston University and the Five College Radio Astronomy Observatory. The survey maps the 
Galactic Ring in the 13 CO J= 1-0 line with an angular and spectral resolution of 46" and 0.2 km s _1 , 
respectively (see Jackson et al. 2006). The observations were performed in both position-switching 
and On-The-Fly mapping modes, achieving an angular sampling of 22". 




Fig. 3. Integrated velocity channel maps of the 13 CO J= 1—0 emission (in green) every ~ 0.6 km s The 
contour levels are 0.3 and 0.8 K km s -1 . Red is the 24 yum emission. 

In our search for molecular gas associated with G26, we analyzed the whole 13 CO J=l-0 data 
cube and found some interesting molecular structures likely related to G26 between 75 and 80 
km s _1 (all velocities are given with respect to the local standard of rest). Figure 3 displays, over 
the MIPS 24 /mi emission (in red), the integrated velocity channel maps of the 13 CO J=l-0 emis- 
sion every ~ 0.6 km s _1 (green with yellow contours), showing the kinematical and morphological 
structure of the molecular gas probably related to the G26 infrared shell. From the panels at 75.0 
and 75.6 km s a molecular clump can be seen, whose peak coincides with the northern brightest 
portion of the infrared shell, where it presents an oblate shape, suggesting an interaction between 
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them. This molecular clump is seen up to ~ 79 km s _1 . In the following panels the molecular gas 
appears to border the northeastern and eastern border of the 24 /im emission, forming an incom- 
plete molecular shell with the lower level contour emission over the IR shell, likely produced by 
the effect of the strong stellar winds, as was found in the LBV star G24.73+0.69 (Petriella et al. 
2012). To better appreciate the molecular gas bordering the north and northeastern portion of the 
24 yum emission, we present in Fig. 4 the 13 CO J=l-0 emission integrated between 77 and 79 km 
s . Additionally, we show in Fig. 5 the 13 CO J=l-0 spectrum obtained toward the peak of the 
molecular feature that it is probably in contact with the northern brightest border of the infrared 
shell (see panels at 75.0 and 75.6 km s _1 in Fig. 3). Clearly, the spectrum is not symmetric and it 
presents a slight spectral shoulder or a less intense component at "redshifted" velocities. It could 
be evidence of turbulent motion in the gas, maybe produced by the strong winds of G26 (see e.g. 
Falgarone et al. 1994). 
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Fig. 5. 13 CO J= 1—0 spectrum obtained toward the peak of the molecular feature that it is probably in contact 
with the northern brightest border of the infrared shell. 
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Taking into account the molecular study presented above, we suggest that the molecular gas is 
very likely related to G26, and hence we adopt the systemic velocity of 77.5 km s for the source, 
deriving a distance of about 4.8 kpc. Hereafter we use this distance for G26. 

From the HI and 13 CO J=l-0 spectra obtained towardxs the source we calculate the total hy- 
drogen column density through Nh = N(HI) + 2N(H2), which is an important factor to compare 
with that obtained from the X-ray emission modeling (see Sect. 3.4). From the HI data we obtain 
N(HI) ~ 1.3 xlO 22 cm 2 . On the other hand, by assuming that the 13 CO J=l-0 line is optically thin 
and using the typical LTE formula as used in Petriella et al. (2012) with T ex = 10 K, we obtain 
N( 13 CO) ~6x 10 15 cm" 2 . Using the relation N(H 2 )/N( 13 CO) ~ 5 x 10 5 (Simon et al. 2001), we 
obtain N(H 2 ) ~ 3 x 10 21 crrT 2 , yielding a total hydrogen column density of Nh ~ 2 x 10 22 cm -2 . 

To roughly estimate some physical parameters of the molecular clump associated with the 
northern maximum at 24 //m of the IR shell (see panels at 75.0 and 75.6 km s _1 ), we again as- 
sumed LTE and T ex = 10 K to derive a column density of N(H 2 ) ~ 1.6 X 10 21 crrT 2 toward this 
structure. The integration was made between 75 and 80 km s . The molecular mass was esti- 
mated by performing the summation of this integrated emission in an elliptical area, centered at / = 
26°479, b = 0°036 with major and minor axes of 55" and 40", respectively and an orientation angle 
of 0°. We obtain a mass of M ~ 125 M , and assuming an ellipsoidal volume we derive a density of 
ne 2 ~ 600 cm 4 , which suggests that the expansion of the infrared nebula is indeed encountering a 
relatively dense molecular clump toward the north. 



3. Multiwavelength analysis towards G26 

The LBV star candidate G26 brights in several wavelengths, which allows us to perform a complete 
study of this object. Figure 6 shows G26 at 24 fim (red) with the radio continuum emission at 20 
cm (in blue with white contours) and some X-ray contours shown in yellow, suggesting that all 
the emissions are related to the same object. Bellow we study the radio continuum, the millimeter 
continuum, and the X-ray emissions in sections 3.1, 3.3, and 3.4, respectively. 




O 
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Galactic Longitude 

Fig. 6. Two-color composite image where the Spfcer-MIPSGAL emission at 24 jim is displayed in red and in 
blue appears the radio continuum emission at 20 cm with white contours with levels of 2.1, 4.1, 6.7, and 9.0 
mjy beairT 1 . The yellow contours show qualitatively the X-ray emission obtained from Chandra observations 
(ObsID = 7493; see Sect. 3.4). The beam of the radio continuum emission appears in the bottom right corner. 
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3.1. Radio continuum emission 

The radio continuum data toward G26 at 1.4 and 4.8 GHz were extracted from the new GPS of 
the Multi- Array Galactic Plane Imaging Survey (Helfand et al. 2006). The emission at 1.4 GHz 
is presented in Fig. 6 (described above). Since the image at 4.8 GHz does not show significant 
morphological difference with that at 1.4 GHz, we do not show it here. 

The dominant feature of the overall morphology is the strong contrast between the brighter 
(W) and much fainter (E) halves of the radio emission. The outermost boundary traces an almost 
elliptical shape with semi-axes of 24" x 18". Most of the radio emission is concentrated in an 
elongated band that runs approximately from north to southeast, which has two radio maximums. 
These maximums are cataloged as the radio sources GPRS 26.470+0.021 and GPRS 26.470+0.025 
(Becker et al. 1994), Rl and R2 in Fig. 6, respectively. The brightest one (Rl), and the only resolved 
with this data set, is coincident with the X-ray emission. It is centered at I = 26°469, b = 0°021 
(18 h 39 m 32 s .2, -05°44'20".7, J2000). The flux densities of Rl at both frequencies were derived by 
integrating the radio emission over a circular area of 11" in diameter defined from the 4.8 GHz 
image, in which the emission is much brighter. We obtained an integrated flux density of 0.025 Jy 
at 1.4 GHz and 0.05 Jy at 4.8 GHz. The other radio source (R2), a point-like one, is located on 
the northern border of the radio emission, at / = 26°470, b = 0°024 (18 h 39 ra 31 s .3, -05°44'12".l, 
J2000). Its peak has a flux density of 0.011 Jy at 1.4 GHz and 0.009 Jy at 4.8 GHz. The errors 
in fluxes are about 20% for both sources at each frequency. In both cases the contributions of the 
extended surrounded emission were subtracted. Thus, we estimate the radio spectral index a (S 
oc y") between both frequencies, obtaining a ~ +0.57 for Rl, which is an expected value for 
thermal emission from the star, and a ~ -0.26 for R2, suggesting a non-thermal radio source. 
On the other hand, to obtain the flux density of the extended radio emission, we subtracted the 
contribution of the two radio sources. We obtain a value of ~0.12 Jy and 0.04 Jy at 1.4 and 4.8 
GHz, respectively. In this case the error in the estimate of the flux densities is about 25% and the 
intrinsic noise of each image was taken into account as well as the uncertainty in the choice of the 

integration boundaries. The radio spectral index of the extended emission turns out to be 0.9, 

clearly non-thermal in nature. 

The observed non-thermal radio continuum should be due to synchrotron emission of relativis- 
tic electrons. The acceleration of the free electrons can be attributed to first-order Fermi accelera- 
tion in shocks within the stellar winds. These shocks can arise either from wind instabilities, or, in 
the case of massive binary systems, close to the contact discontinuity where the stellar winds of two 
stars collide (e.g. Dougherty & Williams 2000 and references therein). It is known that Wolf-Rayet 
stars, which lie on the same evolutionary path as the LBVs, frequently have binary companions, 
which suggests that it would be common that LBV stars form binary systems (Duncan & White 
2002). One possibility is that G26 is a binary system composed of Rl and R2, which are separated 
by ~15" (~0.3 pc or 62 xlO 3 A.U., assuming the distance of 4.8 kpc), and the non-thermal radio 
continuum emission arises from colliding winds. If this is the case, we are indeed looking at a very 
wide binary system (see e.g. Shaya & Oiling 201 1 ; Caballero 2009). Another scenario could be that 
R2 is not a companion star and/or it is just a clump in the non-thermal radio continuum emission. 
In this case, the extended radio emission would be due to wind instabilities from Rl, or colliding 
winds from Rl and an unseen stellar component. 
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3.2. Near-infrared sources 




Galactic Longitude 



Fig. 7. 2MASS JHK three-color image toward G26 (/ in red, H in green, and K s in blue). The white contours 
correspond to the radio continuum emission at 20 cm, whose levels are 2.1, 4.1, and 6.1 mjy beam -1 . 
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Fig. 8. Left: Color-color diagram showing the position of IRSl and IRS2. For comparison we include some 
confirmed LBV stars. The two solid curves represent the location of the main sequence (thin line) and the 
giant stars (thicker line) derived from Bessell & Brett (1988). The parallel dashed lines are reddening vectors 
with the crosses placed at intervals corresponding to five magnitudes of visual extinction. We assumed the 
interstellar reddening law of Rieke & Lebofsky (1985). Right: Color-magnitude diagram of IRSl and IRS2 
The solid curve represents the position of the main sequence at a distance of 4.8 kpc. The reddening vector 
for an 03-type star, with the squares placed at intervals corresponding to five magnitudes of visual extinction, 
is shown with a dashed line. 

To look for a possible stellar companion of Rl, we inspected the Two Micron All-Sky Point 
Source Catalog (2MASS, Cutri et al. 2003) and extracted an image from the 2MASS Image 
Inventory Service. Figure 7 shows a near-infrared (NIR) JHK three-color image toward G26 
with the contours of the radio continuum emission overlaid. We show two infrared sources, 
IRSl (2MASS J 18393224-0544204) coinciding with the radio source Rl, and IRS2 (2MASS 
J18393228-0544148), which is located ~ 5" northeastern from IRSl. In Table 1 we present the 
2MASS photometric data of IRSl and IRS2. The errors in the magnitudes are included between 
parentheses. Both sources appear, in a typical (J-H) versus (H-K s ) color-color diagram (Fig. 8 left), 
as main-sequence stars with a similar visual absorption, suggesting that they could be located at 
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the same distance. For comparison, we include some confirmed LBV stars in this diagram, such as 
HR Car, AFGL 2298, and the Pistol Star. Then, by assuming a distance of 4.8 kpc, we construct a 
K s versus H-K s color-magnitude diagram (Fig. 8 right), which shows that both sources are located 
considerably above the reddening track of an 03V-type star, remarking that they are giant stars, 
mainly IRS1. In a recent study of a stellar cluster located in the W33 complex, Messineo et al. 
(201 1) presented a similar color-magnitude diagram in which blue and red super giants, WR and a 
cLBV star appear quite above the reddening track of an 03V-type star, as in our case. This rough 
NIR photometric study suggests that IRS1 and IRS2 are giant stars probably located at the same 
distance, and hence they could be companions forming a binary system. However, we cannot dis- 
card the hypothesis presented above of the very wide binary system composed of the radio sources 
Rl and R2. The R2 source is associated with the IR source 2MASS J18393143-0544132, which is 
detected only in the K s band (K s = 13.079), and it is impossible to perform photometry. 

Table 1. 2MASS photometric data of IRS 1 and IRS2. 



Source 


J 


H 


K s 


Quality 




(mag) 


(mag) 


(mag) 




IRS1 


7.997(0.021) 


6.526(0.026) 


5.608(0.017) 


AAA 


IRS2 


9.587(0.126) 


8.317(0.100) 


7.572(0.110) 


BAB 



3.3. Millimeter continuum emission 

Millimeter continuum emission from stars may result from thermal free-free, non-thermal syn- 
chrotron, and thermal dust emission (Pallavicini & White 1996). We used the Bolocam Galactic 
Plane Survey (BGPS) to study continuum emission at 1.1 mm toward G26. The BGPS is a 1.1 
mm continuum survey of the Galactic Plane made using Bolocam on the Caltech Submillimeter 
Observatory with 33" FWHM effective resolution (Aguirre et al. 2011). We found a source from 
the BGPS located right upon the LBV star candidate, namely G026.469+00.021, which is dis- 
played in Fig. 9. The positional coincidence between G26 and the millimeter source points to a 
physical connection between them. 

We first explore the possibility that the millimeter continuum emission originates only in the 
dust. Following Rosolowsky et al. (2010), we derive the mass of the emitting dust using 

M = 0A3DF *-^mzi M ®> (i) 

where D is the distance in kpc, F v is the flux at 1 . 1 mm in Jy, and T c { is the dust temperature in K. 
This equation assumes an excitation temperature of 20 K. Taking = 20 K and a distance of 4.8 
kpc, and using the total integrated flux of the source in the 1 . 1 mm band (F v = 204 mJy), we obtain 
a mass of dust of 0.6 M . It is unlikely that this dust component is linked to an ejection event of 
G26 for two reasons. First, as noted by Boyer et al. (2010), cold dust (T<100 K) around LBVs may 
correspond to pre-existing ISM dust swept up by stellar winds. Second, the mass of dust is more 
than an order of magnitude larger than the mass estimated by Clark et al. (2003) (~0.019 M ). 
Moreover, assuming a gas-to-dust ratio of 100, we obtain a gas mass of 60 M , which is greater 
than the mass of the ejection nebulae around other LBVs (see Smith & Owocki 2006 and Clark 
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Fig. 9. Two-color composite image where the Spitzer-MIPSGAL emission at 24 m is displayed in red and the 
smoothed BGPS 1 . 1 mm continuum emission is presented in green with white contours with levels of 100 and 
170 mjy beam -1 . 

et al. 2009). However, in a submillimeter study of r\ Carinae, Gomez (Nee Morgan) et al. (2006) 
estimated a mass of dust between 0.3 and 0.7 M , which may indicate that the gas-to-dust ratio 
in massive stars may differ from the canonical Galactic value. We conclude that if the millimeter 
continuum emission originates in the dust, it is probably be pre-existing ISM dust swept up by 
stellar winds. 

As we showed in Sect. 3.1, the position of G26 coincides with the radio source Rl, which has a 
positive spectral index. This points to a thermal origin for the radio continuum emission. Therefore, 
we investigated the possibility that the millimeter emission originates in the ionized stellar wind 
by means of the thermal free-free mechanism. Then, we analyzed the spectral energy distribu- 
tion (SED) in the radio and millimeter regimes, which relates the flux F v to the wind parameters. 
Assuming that G26 is located at 4.8 kpc, we scaled the stellar parameters derived by Clark et al. 
(2003) using the proportionalities introduced by Hillier et al. (1998): L cc d 2 , M cc d 15 , r eff oc d°, 
where L is the luminosity, M is the mass-loss rate, and T e s is the effective temperature. For the 
luminosity, Clark et al. (2003) quoted log(L/L G ) < 6.0 for a distance of 6.5 kpc, and this value 
drops to ~ 5.4 for a distance of 4.8 kpc. This makes G26 one of the faintest LBVs (see the HR 
diagram of Clark et al. 2009) and a post-red-supergiant-phase candidate (de Jager 1998). For the 
mass-loss rate, the value obtained by Clark et al. (2003), i.e. M = 9.5 x 10~ 5 M /yr, drops to 
5.7 x 10~ 5 M /yr. Finally, the temperature does not scale with the distance. 

After scaling the stellar parameters according to the new distance for G26, we use the following 
equation (Blomme 201 1) to analyze the SED: 

4/3 

F v [mJy] = 0.023-^ (— ) (gyZ 2 v[GHz]f , (2) 

where D is the distance in kpc, M is the mass-loss rate in units of 10~ 6 M Q yr, yt is the mean atomic 
mass, v w is the terminal wind velocity in units of 1000 km s _1 , y is the electron-to-ion density 
number, Z = with aee the helium number abundance, and g = 9.77 [l + 0.13 log 1() (37)] is 

the Gaunt factor, with T the wind temperature in K. The frequency dependence of the Gaunt factor 
is weak and the free-free emission from an ionized stellar wind can be approximated by F v cc y° 6 . 
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v(Hz) 

Fig. 10. Solid lines are the SED of an ionized stellar wind assuming different wind velocities (in km s -1 ). The 
wind parameters used for the calculation are discussed in the text. The dots are the 20 cm (1.4 GHz), 6 cm 
(4.8 GHz), and 1.1 mm (268 GHz) fluxes presented with their errors. 

Assuming a helium abundance = 0.1, we obtain// = 1.3andZ = 1.3. We take the mass-loss 
rate from Clark et al. (2003), scaled to a distance of 4.8 kpc, i.e. M = 5.7 X 10 -5 M Q /yr, as shown 
above. The wind temperature is set equal to the stellar temperature (~ 17,000 K). For a fully ionized 
H and He wind, we obtain y = 1.1. In Fig. 10 we plot the SEDs obtained from eq. 2 for different 
terminal wind velocities (solid lines) and the fluxes in the radio bands at 20 cm (1.4 GHz) and 6 cm 
(4.8 GHz), and the 1.1 mm band (268 GHz). Taking into account that the BGPS has lower angular 
resolution than the radio observations, the 1.1 mm flux may include contributions from both Rl 
and R2 radio sources. However, the expected flux for R2 at 1.1 mm should be ~5 mJy (estimated 
from the spectral index of R2, see Sect. 3.1), which is negligible compared to the total flux (~240 
mJy). Accordingly, we assume that the millimeter emission originates in Rl. From Fig. 10, we see 
that the fluxes from G26 depart from the F v oc v° 6 dependence. Considering the fluxes in the three 
bands, we obtain a ~ 0.37. Only the spectral index between 20 cm and 6 cm (~0.57, see Sect. 3.1) 
approaches the expected value. However, we cannot discard a free-free origin for the millimeter 
continuum emission. 

From this analysis, we conclude that both thermal free-free and dust emission (or a combination 
of them) are plausible mechanisms to explain the presence of BGPS G026. 469+00. 021 toward G26. 

3.4. X-ray emission 

To study the physical connection of the X-ray emission detected from G26 with the observations 
obtained at other wavelengths presented above, we have re-analyzed Chandra observations which 
were first presented by Naze et al. (2012). Our X-ray study improves, with a better fit, the study 
performed by these authors for this particular object. 

3.4.1. X-ray data 

The field of G26 was observed with the ACIS camera of Chandra on 2008 March 10 for 19.5 ks 
(ObsID = 7493). Chandra observations were calibrated using CIAO (version 4.1.2) and CALDB 
(version 3.2.2). To exclude strong background flares that eventually affect the observations, we 
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extracted light curves of photons above 10 keV from the entire field-of-view of the camera, and 
excluded time intervals up to 3<x to produce a GTI file. 

We excluded bad pixels from the analysis using the customary bad-pixel file provided by the 
Chandra X-ray Center for this particular observation. We have searched the data for background 
flares, which are known to affect Chandra data, by examining the lightcurve of the total count rate. 
An additional analysis was performed using the FTOOLS tasks. The spectra were analyzed and 
fitted within XSPECvl 1.0.1 (Arnaud 1996). 

3.4.2. X-ray images 

The sensitivity of the Chandra observations allowed us to detect a medium/hard X-ray point-like 
source at the geometrical center of the G26 infrared nebula. In the soft energy range (i.e. < 1 .2 keV) 
and above the 5.0 keV the source is not detected. This source is located at / = 26?469, b =0°020 
(18 h 39 m 32 s .219, -05°44'19". 15, J2000), and has a signal-to-noise ratio of ~ 7. In Fig. 1 1, we show 
an ACIS image in the energy band of 1-5.5 keV with radio contours at 1.4 GHz superimposed. As 
can be seen, the Rl peak is coincident with the point X-ray source detected by Chandra, and no 
X-ray emission is detected at the position of R2. To extract counts in the source region and compute 
net counts we used the dmextract tool. As a result, 400 counts were computed. 




Galactic Longitude 



Fig. 11. Chandra ACIS image, with a size of 0.4 x 0.4 arcmin, of G26 in the 1.0-5.0 keV energy band. The 
smoothed image was convolved with a Gaussian function of a kernel radio of 3x3 pixels (~ 1.5 arcsec). The 
green contours represent the radio continuum emission at 20 cm with levels of 2.1, 4.1, and 6.1 mjy beam -1 . 

Finally, to search for variability in the ACIS-I observation, we used the photon arrival times in 
the 1.0-5.0 keV band and the CIAO "glvary" tool using the Gregory-Loredo algorithm (Gregory & 
Loredo 1992). No hints of variability were detected during this observation. 

3.4.3. Spectral X-ray analysis 

We extracted the X-ray spectrum of G26 from a circular region with a radius of 4 arcsec using the 
CIAO tool specextract. The background spectrum was estimated from an annular region with radii 
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Energy (keV) 

Fig. 12. Chandra X-ray spectrum of G26 in the 1.0-5.5 keV. The solid lines indicate the best-fit (VAPEC) 
model (see Table 2). Lower panel: Chi squared residual of the best-fit model. 



of 4 and 10 arcsec. The spectrum is grouped with a minimum of 16 counts bin 1 and the^ 2 statistic 
is used. The errors quoted are 90%. 

The spectrum is shown in Fig. 12 and it is well represented by a variable APEC model 
(VAPEC), which provides the best acceptable fit. The X-ray parameters for the best fit are given in 
Table 2. The spectrum is dominated by atomic emission lines of Si and S, a temperature of ~ 0.89 
keV, and a high neutral hydrogen absorption column. These characteristics suggest that the X-ray 
emission has an optically thin thermal plasma origin. Unfortunately, the poor quality of the data 
above 4 keV does not allow us to check for the FeXXV (6.4-6.7 keV) complex lines that could 
help us to constrain the thermal and/or non-thermal emission at the hard part of the spectrum. For 
references about the detection of the FeXXV complex lines in massive stars see e.g. Hyodo et al. 
(2008), Albacete Colombo et al. (2007), and Smith et al. (2004). The obtained absorbing column 
density agrees well with that obtained from the CO and HI data (see Sect. 2.2), confirming therefore 
the amount of material between us and the source. This fact supports the obtained radial velocity 
for the source and indirectly confirms the derived distance. 

We calculated the unabsorbed X-ray luminosities in the 1.5-5 keV band for a distance of 4.8 
kpc, obtaining Lx ~ 4.7 x 10 33 erg s" 1 , and yielding log[Lx/LBoiJ = -5.35 + 0.05. As mention 
above, the spectral X-ray analysis reveals that the emission has an optically thin thermal origin, 
which is dominant in the medium and hard energy range. These characteristics seem to be indicat- 
ing wind-wind collision shocks from a binary system of massive stars (e.g. De Becker et al. 2005; 
Pittard & Parkin 2010), reinforcing the hypothesis presented in Sees. 3.1 and 3.2. This result does 
not seem to be common in LBVs. At present only very few X-ray emissions toward LBVs indicate 
wind-wind collision shocks (see Naze et al. 2012). Indeed, as these authors point out, the LBVs 
as a class are clearly not bright X-ray emitters, and the few cases of X-ray detections should have 
an extrinsic cause, such as binarity. If it is indeed common that LBVs possess companions like 
WRs do, following Naze et al. (2012), there are basically two possibilities to explain the observed 
X-ray emission toward these objects: (i) the companion is quite close to the LBV, or (ii) the com- 
panion is relatively distant from the LBV. In the first case, the intrinsic emission of the companion 
may be hidden by the strong absorption of the dense wind, but then the conditions are favorable 
for an X-ray bright wind-wind collision. In the second case, no emission from a wind-wind colli- 
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Table 2. X-ray spectral parameters for G26 



Model 




Parameters 




PHABS 




N H [cm- 2 ] 


3.46(±0.43)xl0 22 


VAPEC 


kT [keV] 


0.89± 0.12 


[Si/H] 


0.68± 0.21 


[S/H] 


0.97+ 0.42 


Norm 


2.04(±1.6)xl0- 3 


^/d.o.f. 


0.86/55 


Flux(1.0-2.5)[ergcrrr 2 s" 1 ] 


1.69(±0.02)xl0" 12 


Flux(2.5-5.0)[erg cm" 2 s _1 ] 


1.02(±0.03)xl0- 13 


Total Flux( 1.0-5.0) [erg crrT 2 s" 1 ] 


1.79 (±0.06)xl0~ 12 



Normalization is defined as 10~ 14 /47rD 2 x J"n w n e dV, where D is distance in [cm], n H is the hydrogen density 
[cm ~ 3 ], n e is the electron density [cm -3 ], and V is the volume [cm 3 ]. The flux in the two energy ranges is 
absorption-corrected. Values in parentheses are the single parameter 90% confidence interval. The abundance 
parameter is given relative to the solar values of Anders & Grevesse (1989). 

sion is expected due to wind dilution, but the intrinsic emission of the companion would be easily 
detectable, because the tenuous wind of the LBV cannot hide it anymore. High sensitivity obser- 
vations are needed to explore these possibilities and the limit between them. If we consider this 
statement to be correct, the nature of the X-ray emission from G26 would discard the very wide 
binary system formed by Rl and R2 (see Sect. 3.1). 

Finally, it is interesting to mention that G26 was first detected by the ASCA telescope (source 
AX J183931-0544) by Sugizaki et al. (2001). In that work the X-ray flux is a factor 2 greater than 
the flux we obtained. However, this comparison should be taken with caution, no details about the 
model used to compute the ASCA flux and other parameters of the best-fit. Therefore, we cannot 
confirm a variable behavior of G26. 

4. Summary 

The luminous blue variable (LBV) stars are peculiar very massive stars. The study of these stellar 
objects and their surroundings is important for understanding the evolution of massive stars and 
the impact that they produce in the interstellar medium. We have performed a multiwavelenght 
study of the LBV star candidate G26.47+0.02 (G26) and its environment. The main results can be 
summarized as follows: 

(1) G26 presents an infrared nebula that brights at 24 //m. We found that this nebula is partially 
surrounded by a molecular shell (seen in the 13 CO J=l-0 line), suggesting an interaction between 
the strong stellar winds from G26 and the molecular gas. From the HI absorption and the molecular 
gas study we conclude that G26 is located at a distance of ~ 4.8 kpc. 

(2) Scaling the stellar parameters derived in previous works to the distance of 4.8 kpc, we 
found that G26 would be one of the faintest LBVs, suggesting it to be a post-red-supergiant-phase 
candidate. 
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(3) The radio continuum emission at 6 and 20 cm shows two radio sources (Rl and R2) and 
extended emission around them within the G26 infrared nebula. Source Rl coincides with the 
central stellar object G26, and R2 is ~ 15" towards the north. From the radio spectral index study 
we found that Rl presents thermal emission, while R2 and the extended emission are non-thermal. 
We suggest that the non-thermal radio emission might arise from the stellar wind instabilities, or, in 
the case of a massive binary system, from wind-wind collisions. In this scenario, Rl and R2 could 
be companions in a binary system. 

(4) From the Two Micron All-Sky Point Source Catalog we found two IR sources with similar 
photomectrical characteristics, one coinciding with Rl, and the other one ~ 5" toward the north- 
east. A rough JHK photometric study suggests that both sources are giant stars probably located at 
the same distance. We present it as another scenario for a possible binary system. 

(5) We found that the 1.1 mm continuum Bolocam source G026. 469+00. 021 coincides with the 
G26 central source. We conclude that both thermal free-free and dust emission (or a combination 
of them) are plausible mechanisms to explain the millimeter continuum emission. 

(6) Chandra observations reveal a medium/hard X-ray point-like source at the geometrical cen- 
ter of G26, coinciding with the radio source Rl . No X-ray emission is detected at the position of the 
radio source R2. No hints of variability were detected during the Chandra observation. Our results 
in the X-ray regime improve, through a better fit, the results obtained by Naze et al. (2012). The 
spectral X-ray analysis reveals that the emission has an optically thin thermal origin, dominant in 
the medium and hard energy range. These characteristics seem to be indicating wind-wind collision 
shocks from a binary system of massive stars. 

Our main result points to the hypothesis that G26 could be a binary system. We give two 
possible scenarios for that: (i) a very wide binary system composed of the radio sources Rl and 
R2, and (ii) a binary system composed of the infrared sources IRS1 (coniciding with Rl) and IRS2. 
With the present data we cannot exclude either of them. Binary systems appear to be common in 
WR stars and thus probably in LB Vs. Therefore multiwavelength and long-term observations are 
needed to detect some possible variable behavior to confirm the binary nature of the system and to 
reliably determine its components. 
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